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1. Introduction. 

For a finite-dimensional integrable Hamiltonian systems (FDIHS), let m denote the 
number of degrees of freedom, and Pi,i = l,...,m, be functionally independent inte- 
grals of motion in involution, the separation of variables means to construct m pairs of 
canonical separated variables Vk,Uk,k = 1, m, [1,2,3] 

{uk,ui} = {vk,vi} = 0, {vk,ui} = Ski, k,l = l,...,m, (1.1) 

and m functions fk such that 

fkiuk,Vk,Pi,...,Pm) ^Q, /c = l,...,m, (1.2) 

which are called separated equations. The equations (1.2) give rise to an explicit fac- 
torization of the Liouville tori. For the FDIHSs with the Lax matrices admitting the 
r-matrices of the XXX, XX Z and XYZ type, there is a general approach to their sep- 
aration of variables [1-6]. The corresponding separated equations enable us to express 
the generating function of canonical transformation in completely separated form as an 
abelian integral on the associated invariant spectral curve. The resulted linearizing map 
is essentially the Abel map to the Jacobi variety of the spectral curve, thus providing a 
link with the algebro-geometric linearization methods given by [7-9]. 

The separation of variables for a FDIHS requires that the number of canonical sep- 
arated variables Uk should be equal to the number m of degrees of freedom. In some 
cases, the number of Uk resulted by the normal method may be less than m and one 
needs to introduce some additional canonical separated variables. So far very few mod- 
els in these cases have been studied. These cases remain to be a challenging problem 
[3]. 

The separation of variables for constrained flows of soliton equations has been stud- 
ied (see, for example, [4,10-14]). In recent years binary constrained flows of soliton 
hierarchies have attracted attention (see, for example, [15-22]). However the separation 

of variables for binary constrained flows has not been studied. The degree of freedom 

for high-order binary constrained flows admitting 2x2 Lax matrices M = 

is a natural number 2N + ko. Via the Lax matrix M, N + ko pairs of canonical separated 
variables ui, ...,UN+ko be introduced by the set of zeros of -B(A) and Vk — A{uk), 
and N + ko separated equations can be found from the generation function of intrgrals 
of motions. In previous papers [23,24] we pesented a method with two different ways for 
determining additional N pairs of canonical separated variables and additional N sepa- 
rated equations for flrst binary constrained flows with 2N degree of freedom. The main 
idea in [23,24] is to construct two functions -B(A) and ^(A) and deflne un^i, ...,U2n 
by the set of zeros of B{X) and VN^k — A{uN^k)- The ways for constructing B{X) 
and ^(A) in [23] and [24] are some different. In present paper we propose a completely 
different method from that in [23,24] to introduce the additional N separated variables 
and N separated equations for high-order binary constrained flows with 2N + ko de- 
gree of freedom. It is observed that the introduction of Vk has some link with integrals 



A{X) B{X) \ 
IC(A) -AiX) 
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of motion and should lead to the separated equations. We find that there are N in- 
tegrals of motion Pjv+fco+i) ■■■tP^N+ko among the 2N + ko integrals of motion for the 
high-order binary constrained flows which commute with A{X) and B{X). This ob- 
servation and property stimulate us to directly use the additional integrals of motion 
to define both the pairs of additional separated variables and N separated equa- 
tions by VN+ko+j = PN+ko+j,j — 1,...,N. Then we can find the conjugated variables 
UN+ko+ji 1) ■■■) ^1 commuting with A{X) and B{X). In contrast to the method in [23,24], 
this method is easer to be applied to the high-order binary constrained flows. 

We will also present the separation of variables of soliton equations. The first step is to 
factorize (1 -|- 1)— dimensional soliton equations into two commuting x— and FDIHSs 
via high-order binary constrained fiows, namely the x— and t— dependences of the soliton 
equations are separated by the x— and t— FDIHSs obtained from the x- and t-binary 
constrained fiows. The second step is to produce separation of variables for the x— and 
t— FDIHDs. Finally, combining the factorization of soliton equations with the Jacobi 
inversion problems for x— and t— FDIHSs enables us to establish the Jacobi inversion 
problems for soliton equations. If the Jacobi inversion problem can be solved by the 
Jacobi inversion technique [7], one can obtain the solution in terms of the Riemann- 
theta function for soliton equations. We illustrate the method by KdV, AKNS and 
Kaup-Newell (KN) hierarchies. The paper is organized as follows. 

In section 2, we first recall the high-order binary constrained fiows and factorization of 
the KdV hierarchy. Then propose the mehtod for introducing the N pairs of additional 
separated variables. We illustrate the method by both first binary constrained fiow 
and second binary constrained fiow. Finally we present the separation of variables for 
KdV hierarchy. In section 3 and 4, the method is applied to the AKNS hierarchy and 
KN hierarchy, respectively. In fact this method can be applied to all high-order binary 
constrained fiows and other soliton hierarchies admitting 2x2 Lax pairs. 

2. Separation of variables for the KdV equations. 

We first recall the high-order binary constrained fiows of the KdV hierarchy. 

2.1 High-order binary constrained flows of the KdV hierarchy. 

Consider the Schrodinger equation [25] 



(t>xx + {X + u)(j) = 



which is equivalent to the following spectral problem 




(2.1) 



Take the time evolution law of as 



0,„=W-)(k,A)0, 



(2.2) 
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where 



ao = 60 = 0, Co = -1, ai = 0, 61 = 1, 

bk+i = Lbk = --L''~^u, , ak = --bk,x, Ck = —-bk,xx - h+i - hu, k = l,2,---, 

L = —d'^ -u+^d-^Ux, d = dx, d-^d = dd'^ = 1. (2.3) 

The compatibihty condition of (2.1) and (2.2) gives rise to the n-th KdV equation 
which can be written as an infinite-dimensional Hamiltonian system [25] 

ut^ = -2bn+2,x = dL^u = (2.4) 

with the Hamiltonian ff„ = and ^ = -2&„+2. 

For 77. = 1 we have 



0,, = V^W(«,A)0, F« = ( 



\u\+\uxx+\u'^ -\ux 



(2.5) 



and the equation (2.4) for n = 1 is the well-known KdV equation 



ut^ = -^{uxxx + Quux). (2.6) 



The adjoint spectral problem reads 



By means of the formula in [26] , we have 

'(piijji (I)iip2 \ dU{u,X) 



According to [15-22], the high-order binary x-constrained flows of the KdV hierar- 
chy (2.4) consist of the equations obtained from the spectral problem (2.1) and the 
adjoint spectral problem (2.7) for N distinct real numbers Xj and the restriction of the 
variational 

derivatives for the conserved quantities H^^ (for any fixed ko) and A^-: 

*l,a; = ^2, ^2,x = -A$i - u^i, (2.8a) 
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= A*2 + li*2, *2,x = -*l, (2.8b) 

^-E^ = -2W<*2,$i>=0. (2.8c) 

Hereafter we denote the inner product in by < ., . > and 

= (0ii, • • • 50iiv)^, *i = (V'ii) • • • ,'0iAr)^, *=1,2, A = diag{Xi,- ■ ■ ,Xn)- 

The binary t^-constrained flows of the KdV hierarchy (2.4) are defined by the replicas 
of (2.2) and its adjoint system for N distinct real number Xj 

(2.9a) 

as well as the n-th KdV equation itself (2.4) in the case of the higher-order constraint 
for ko > 1 

Ut^ = -2bn+2,x. (2.9b) 

(1) For ko = 0, we have 

&2 = -^ii= ^ < *2,$i >, i-e., w = -<*2,$i>. (2.10) 

By substituting (2.10), (2.8a) and (2.8b) becomes a finite- dimensional Hamiltonian sys- 
tem (FDHS) [18] 

dFi ^ dFi ^ dFi ^ dFi 

$1. = ^, ^2. = ^, = = ^2.11) 

i^l =< *1, $2 > - < A*2, $1 > + J < *2, $1 >^ • 

Under the constraint (2.10) and the a;-FDHS (2.11), the binary ti-constrained flow 
obtained from (2.9a) with V^'^^ given by (2.5) can also be written as a ti-FDHS 

5F9 5F9 dF-y 

5*2' 9^2' ^ ^ 

F2 = - < A^^'2, $1 > + < A^'l, $2 > +^ < *2, $1 >< A^'2, $1 > 

1 1 

+ - < *2,$1 >< *1,$2 > +„(< *2,$2 > - < >)^. 

2 o 
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The Lax representation for the a;-constrained flow (2.8) and the t„-constrained flow 
(2.9) can be deduced from the adjoint representation of (2.1) and (2.2) by using the 
method in [27,28] 

M, = [C/,M], Mt^ = [V^^\M], (2.13) 

where U and F^") are obtained from U and V^'^^ under the system (2.8), and the Lax 
matrix M is of the form 

/A(A) 5(A) \ 
\C{X) ~A{X))- 

The Lax matrix M for x-FDHS (2.11) and ti-FDHS (2.12) is given by 

C(A) = -A + l<<l?..*,>+iX:t^- (2-") 

The generating function of integrals of motion for (2.11) and (2.12) yields 

^^(A) + B{\)C{X) ^ P(A) = -A + f^ij^. + ^X-Xr ^' ^^-^^^ 
where Pi, ...,P2n are independent integrals of motion for the FDHSs (2.11) and (2.12) 

Pj = \M23 + {-\Xj + \< *2,$1 >)^2i01,- 

+ 8 XI A - - Afe ~ V'2j02j)(V'ifc0ifc - '4^2k<p2k) + 4V'ii02jV'2fc0ifc], (2.16a) 



PN+j = ^{'(piAij + V'2j</'2i), j = 1, N. (2.16b) 



It is easy to verify that 

JV N 



F, = 2j2Pj, F2 = 2j](A,P,+P^+,.)- (2.17) 

i=i i=i 



With respect to the standard Poisson bracket 

JV 



fj g^ _ ^9 ^ df dg df dg df dg 
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by calculating the formulas like (2.31), it is easy to verify that 

{A\X) + B{\)C{\),A\fi) + B{fi)C{ii)} = 0, (2.19) 

which implies that Pi, .., P2N are in involution, (2.11) and (2.12) are FDIHSs and com- 
mute with each other. The construction of (2.11) and (2.12) ensures that if (^^i, ^^2, ^i, ^2) 
satisfies the FDIHSs (2.11) and (2.12) simultaneously, then u defined by (2.10) solves 
the KdV equation (2.6). 
Set 

A\X)+B{X)C{X) = Xj2Fk^~\ (2-20) 

A;=0 

where F^, k = 1,2, are also integrals of motion for both the x-FDHSs (2.11) and the 
t^-binary constrained flows (2.9). Comparing (2.20) with (2.15), one gets 

TV 

Fo = -l, Fi=0, Ffe = ^[A)-2p, + (/c-2)A^"-3p2+,.], A; = 2, 3,.... (2.21) 

By employing the method in [28,29], the t„-FDIHS obtained from the t„-binary con- 
strained flow (2.9) is found to be of the form 

dFn+i _ dFn+i _ dFn+i _ dFn+i , , 

m=0 li + ...+lm+i=n+2 

where h > 1, Im+i > 1, cko = I5 = |) 0:2 = §, and [28,29] 

m—2 ^ m— 1 

«m = 2q;^_i + ^ aiam-i-1 ~ 2 ^ oiio^m-l, m>3. (2.22c) 
1^1 1=1 

The n-th KdV equation (2.4) is factorized by the x-FDIHS (2.11) and the tn-FDIHS 
(2.22). 

(2) For ko = 1, one gets 

b3=^{Ua:x + 3u^) = \<'^2,^l> . (2.23) 

By introducing q — u,p — \ux, (2.8a), (2.8b) and (2.23) can be written as a x-FDHS 
dFi dFi dFi dFi 
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Fi = - < A*2, $1 > + < $2 > -g < *2, $1 > +2p^ + ^q^. 
Under the constraint (2.23), V^^^ becomes 

V^'^=( ^2 1 . ^ . 1 2 ^"^'V (2.25) 

Under the constraint (2.23) and the x-FDHS (2.24), the binary ti -constrained flow 
consists of (2.9a) with V^^^ replaced by V^^^ and (2.9b) given by (2.6) can also be 
written as a ti-FDHS 

dF2 dF2 dF2 dF2 

*.,.=-g^, i=1.2. Pn=-^. (2.26) 

F2 = - < A^*2, $1 > + < A^'i, $2 > ~q < A*2, $1 > ~q < $2 > 

1 9 1 9 

+p < *i, $1 > -p < ^'2, $2 > +2 < $1 > -^q < ^'2, $1 > • 

The Lax representations for the x-FDHS (2.24) and the ti-FDHS (2.26), which can 
can be deduced from the adjoint representation of (2.1) and (2.2), are given by (2.13) 
with V^^^ defined by (2.25) and U obtained from U by using q instead of u as well as 
M given by 

3 = 1 ^ 

The generating function of integrals of motion for (2.24) and (2.26) yields 

iV p p2 

A\X) + B{X)C{X) ^ P(A) = -A^ + Po + E[ + (A"=aJf^' ^^'^^^ 

where Pq, ...,P27V are independent integrals of motion for the FDHSs (2.24) and (2.26) 

and Po = l-P^i, 

Pj = -^A|V'2i01i + ^AjVlj02j - ^Aj#2i01i - ^q'iplj(p2j 
+ \p{'^lj4>lj - •02i02i) + ^(< ^'2, $1 > -^q^)'ip2j(j)lj 



A METHOD TO INTRODUCE ADDITIONAL SEPARATED VARIABLES 



9 



+ o 5^ ^[(V'lj^lj - 'ip2j(p2j){lplk(plk - ^2fc02fc) + 4V'li02j^2fc01fc], (2.29a) 

Pn+j = + V'2,</'2i), j = 1, N. (2.29b) 

We have 

N 

Fi = 2Po, F2 = 2^P,. (2.30) 

i=i 

Similarly, it can be shown that (2.24) and (2.26) are FDIHSs and commute with each 
other. The KdV equation (2.6) is factorized by x-FDIHS (2.24) and ti-FDIHS (2.26). If 
(*i,\&2,P, $i,$2,?) satisfies the FDIHSs (2.24) and (2.26) simultaneously, then u = q 
solves the KdV equation (2.6). 

2.2 The separation of variables for the KdV equations. 

(1) For the case kg = 0, we first consider the separation of variables for FDIHSs 
(2.11) and (2.12). With respect to the standard Poisson bracket (2.18), it is found that 
for the A{X) and B{X) given by (2.14) we have 

{A{X),A{t,)} = {B{X),B{,,)} = 0, {A{X),B{,,)} = ^^^^[B{^^)-BiX)]. (2.31) 

An effective way to introduce the separated variables Vk,Uk and to obtain the separated 
equations is to use the Lax matrix M and the generating function of integrals of motion. 
The commutator relations (2.31) and the equation (2.15) enable us to define the first 
N pairs of the canonical variables tti, ...,un by the set of zeros of B{X) [1-3] 

where 

N N 

R{X) ^l[{X-Uk), K{X) = J] (A - Afe), 

k^l 

and vi, vn by 

Vk = 2A{uk), k = l,...,N. (2.32b) 

The commutator relations (2.31) guarantee that ui,...,un and vi,...,vn satisfy the 
canonical conditions (1.1) [1-3]. Then substituting Uk into (2.15) gives rise to the first 
N separated equations 



Vk 



2^(«fc) = 2/PK) = 2, 



. -«^ + E[^l- + r^%^]' k = l,...,N. 

(2.33) 
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The FDIHSs (2.11) and (2.12) have 2N degrees of freedom, we need to introduce the 
other N pairs of canonical variables VkiUkik = N -\- 1, ...^ 2N. Notice that Pn+j given 
by (2.16b) are integrals of motion for the FDIHSs (2.11) and (2.12), and satisfy 

{B{X), Pn+j} = {A{X), Pn+j} = 0. (2.34) 

Thus we may define 

VN+j = 2PN+j = ^(V'ij0ij + V'2j02j), j = 1, N, (2.35a) 
which also give rise to the separated equations. It is easy to see that if we take 

UN+j=ln^, j = l,...,N, (2.35b) 

then 

{vN+j:UN+k} = Sjk, {vN+j,VN+k} = {uN+j,UN+k} = 0, j,k = l,...,N, (2.36) 

{B{X),UN+j} = {A{X), UN+j} = {B{X), VN+j} = {A{X),VN+j} = 0. (2.37) 

We have the following proposition. 

Proposition 1. Assume that Xj , 4>ij , ipij e R, i = l,2,j = 1,...,N. Introduce the 
separated variables ui,...,U2n and vi,...,V2n by (2.32) and (2.35). If ui, ...jUn, are 
single zeros of B{X), then vi,...,V2n and ui,...,U2n are canonically conjugated, i.e., 
they satisfy (l-l). 

Proof. By following the similar method in [1-6,23,24], it is easy to show that t"!, 
and Ml, ...,un satisfy (1.1). Notice B'{u}S) ^ 0. Hereafter the prime denotes the differ- 
entiation with respect to A. It follows from (2.36) and (2.37) that 

= {uN+k-, B{ui)} = B'{ui){uN+k, ui} + {uN+k, B{ii)}\i^=u^ = B'{ui){uN+k, 

{Vk,UN+l} = 2{A{uk),UN+l} 

= 2A'{uk){uk, UN+i} + {A{X), UN+i}\x=Uk = '2A'{uk){uk, un+i}, (2.38) 

which leads to {itjv+fc)ii«} = {uN+ki^i} = 0. Similarly we can show that {vN+k,ui} = 
{vN+k,vi} = 0. These together with (2.36) complete the proof. 
It follows from (2.32a) and (2.35b) that 

AT 

u= -< ^2,^1 >=2^{uj - Xj), (2.39) 
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or 



The separated equations are given by (2.33) and (2.35a). Replacing by the partial 
derivative of the generating function S of the canonical transformation and inter- 
preting the Pi as integration constants, the equations (2.33) and (2.35a) give rise to the 
Hamilton-Jacobi equations which are completely separable and may be integrated to 
give the completely separated solution 

S{m, ...,U2n) = Y1U 2,/P{X)dX + 2PN+kUN+k]- (2.41) 
k=i 

The linearizing coordinates are then 

dS ^ f^^ Pat ■ 

QN+i = = 2 E / dX + 2uN^i, i = 1, N. (2.42b) 

9PN+i {X-Xify/PiX) 

By using (2.17), the linear flow induced by (2.11) is then given by 

dFi dFi 
Qi = li + x—— ^^i + 2x, QN+i = '^lN+i + x— =27Ar+i, z = 1, AT. (2.43) 

Hereafter 7j, i = 1, 2N, are arbitrary constants. Combining the equation (2.42) with 
the equation (2.43) leads to the Jacobi inversion problem for the FDIHS (2.11) 



k=l 

N 



^1 ix-K)vm ''^"^'''- *^ 



fe=i 



Uk p 



dX + UN+i = ^N+i, i = l,...,N. (2.44b) 



(A-A02/P(A) 



If, by using the Jacobi inversion technique [7], ■4>2ji < '^2^^i > given by (2.39) and 
(2.40) can be obtained from (2.44), then (t>2j,'ipij can be found from the first and the 
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last equation in (2.11) by an algebraic calculation, respectively. The i4>ij , 4'2j , ipij , 'il^2j) 
provides the solution to the FDIHS (2.11). 

By using (2.17), the linear flow induced by (2.12) is then given by 

dF2 

Qi = li + -gp:^^ = 7i + 2Ajti, 

dFo 

Qn+i = 2^N+i + 777^ ^1 = 27jv+i + 4Pjv+iti, i = 1, iV, (2.45) 

where 7^ are arbitrary constants. Combining the equation (2.42) with the equation 
(2.45) leads to the Jacobi inversion problem for the FDIHS (2.12) 



fe=l 



/ Pat 
-Ii±l—=d\ + UN+i = iN+i + 'iPN+itu i = 1, N. (2.46b) 

Finally, since the KdV equation (2.6) is factorized by the FDIHSs (2.11) 
and (2.12), combining the equation (2.44) with the equation (2.46) give rise to the 
Jacobi inversion problem for the KdV equation (2.6) 



77 —-j==d\^^i + 2x + 2Xitu i = l,...,iV, (2.47a) 



N 



{X-Xi)^/P{Xj 



E/ 7^ -^^^^J==dX + UN+^=lN+^ + 2PN+^tl, i = 1, TV. (2.47b) 

t^i-i {x-Xiy^jp{X) 

If, by using the Jacobi inversion technique [7], u given by (2.39) can be found in terms 
of Riemann theta functions by solving (2.47), then u provides the solution of the KdV 
equation (2.6). 

In general, since the n-th KdV equation (2.4) is factorized by the x-FDIHS (2.11) and 
the tn-FDIHS (2.22), the above procedure can be applied to find the Jacobi inversion 
problem for the n-th KdV equation (2.4). We have the following proposition. 

Proposition 2. The Jacobi inversion problem for the n-th KdV equation (2.4) is given 
by 



^ r'^k 1 

1 1 *f 



, d\ — 'ji + 2x 

(A-A,)yP(A) 



+*nE(9r"'«- E >^r'-'Fi,...Fi^, i = i,...,N, 



2 

m=0 h + ...+lm+i=n+2 
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iV 

E 



N+i 



\J (A-A,)2/P(A) 



dX + UN+i = iN+i 



m=0 ii + ...+Z^+i=n+2 

where Zi > 1, ...,Zrn+i > 1 ond Fcj, ■■■Fi^, are given by (2.21). 

(2) For the case /cq = 1, we now consider the separation of variables for FDIHSs 
(2.24) and (2.26). With respect to the standard Poisson bracket, it is found that the 
A{X) and -B(A) given by (2.27) also satisfy commutator relation (2.31). In the same 
way, the first + 1 pairs of the canonical variables wi, un+i can be introduced by 
the set of zeros of B{X) 



N 



B{X) = X--q+-J2 



1 1 ^ ^^2,01, _ RiX) 



2^ ' 2^ A-A, KiXy 



(2.48a) 



where 



N+l 



N 



R{x) = n - «'^)' ^(^) = n - ^fc)' 



k=l 



k=l 



and vi, ...,vn+i by 



Vk^2A{uk). k = l,...,N+l. (2.48b) 
Then substituting into (2.28) gives rise to the first A'" + 1 separated equations 



Vk 



2A{uk) = 2^yP{uk) = 2 



i 



N 



-u) 



p2 



^ Wfe - Xj {uk - Xj) 



2J' 



k = l,...,N+l. (2.49) 
The additional pairs of canonical variables can also be defined by the same way 



VN+l+j = 2Pjv+i = ^(V'lj^li + V'2j02j), j = 1, 



(2.50a) 



MAr+i+j = ln— , j = l,...,N. 



(2.50b) 



In the same way we can show the following proposition. 
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Proposition 3. Assume that Xj,4>ij,ipij e R, i = 1,2, j = 1,...,N. Introduce the 
separated variables ui, ...,U2N+i and vi, ...,V2n+i by (2.48) and (2.50). If ui, ...,un+i, 
are single zeros of B{X), then vi, ...,V2n+i and ui, ...,U2n+i o-re canonically conjugated, 
i.e., they satisfy (l-l). 

It follows from (2.48) and (2.50) that 

iV+l N 

u = q = 2^Uj - 2^Xj, (2.51a) 

3=1 j=i 



'2i?(Aj)e"^+i+^- _ /2i?(Aj)e-"^+i+j- 



'^j = \l ^77V-^ ' V's,- = W ^T^TTT^ > J = l,-,^- (2.51b) 



The separated equations (2.49) and (2.50a) may be integrated to give the completely- 
separated solution for the generating function S of the canonical transformation 

iV+l .ufc N 

Siui,...,U2N+i)=J2 2y/P(XjdX + 2j2PN+kUN+i+k, (2.52) 

k=l k=l 

where -P(A) is given by (2.28). 

In the exactly same way, one gets the Jacobi inversion problem for the FDIHS (2.24) 

N+l I 

E / ^H7Tf^^ = ^o + 2x, (2.53a) 
fe=i V^K^) 

N+l Uk I 

Ylij (A - A ^ UN+i+i = iN+i, i = l,...,N, (2.53c) 

the Jacobi inversion problem for the FDIHS (2.26) 

iV+l 



k=l 

N+l 



/ Pn ■ 
—^^y==dX + UN+i+i = ^N+i, i = l,...,N. (2.54c) 
(A-A^)V^(A) 



Finally we have the following proposition. 
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Proposition 4. The Jacobi inversion problem for the KdV equation (2.6) 

N+l 



E 

k=l 



AT+l 

E 

fe=l 



AT+l 

E 

k=l 



Uk 



(A-Ai)A/P(A) 



Uk 2 

-^_dA = 70 + 2x, (2.55a) 
rfA = 7i + 2ii, i = l,...,N, (2.55b) 



N+i 



dX + UN+i+i^^N+i, i = l,...,N. (2.55c) 



If, by using the Jacobi inversion technique [7], u given by (2.51a) can be found in 
terms of Riemann theta functions by solving (2.55), then u provides the solution of the 
KdV equation (2.6). 

In general, since the n-th KdV equation (2.4) is factorized by the x-FDIHS (2.24) 
and the ^n-FDIHS obtained from (2.9) under (2.24), the above procedure can be applied 
to find the Jacobi inversion problem for the n-th KdV equation (2.4). 

(3) The method can be applied to all high-order binary constrained flows 
(2.8) and (2.9) as well as the whole KdV hierarchy. For any fixed ko, 
by introducing the so-called Jacobi-Ostrogradsky coordinates, the high-order binary 
x-constrained fiow (2.8) can be transformed into a x-FDIHS with degree of freedom 2iV-|- 
ko- Under the x-FDIHS, the binary t^-constrained fiow (2.9) can also be transformed 
into a t„-FDIHS. The Lax representation for the x- and t„-FDIHS can be deduced from 
the adjoint representation of (2.1) and (2.2) by using the method in [27,28]. By means 
of the Lax matrix we can introduce the first N + ko canonical variables ui, ...^uiyj^ko 
by the set of zeros of B{X) and — 2A{uk), k = 1, N + ko. Then the additional A'" 
canonical separated variables can be defined by 

VN+ko+j = SPat+j = ^(V'lj^li + V'2j02j), UN+ko+j = j = 1' 

Finally, since the n-th KdV equation (2.4) is factorized by the a;-FDIHS and the 
FDIHS, in the exactly same way we can obtain the Jacobi inversion problem for (2.4). 
The above scheme can be applied to all soliton equations admitting 2x2 Lax pairs. 

3. The separation of variables for the AKNS equations. 

3.1 Binary constrained flows of the AKNS hierarchy. 

For the AKNS spectral problem [30] 
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Take 



where 



(3.2) 



(3.3) 



ao = -1, 6o = Co = 0, ai = 0, 6i = g-, ci = r, a2 = ^qr, 

(fefe+l) ^ ^(fefc)' ak ^ d~^{qck - rbk), k=l,2,---, 

J _ ^ ( d — 2rd~^q 2rd~^r \ 
2 V -'2qd-\ -d+2qd-^r J ' 

The AKNS hierarchy associated with (3.1) and (3.2) reads [30] 

-0 "'-"^V 

For 77. = 2 we have 

and the AKNS equation (3.4) for n = 2 reads 

1 1 
9*2 = -i^Qxx + = -ra,x - r^q. (3.6) 

The adjoint AKNS spectral problem is of the same form as equation (2.7). 
We have 



5u~ ^""^(,02^1 ^^^2 ) du " (,^^2^ J' ^^■'^^ 

which should be read componentwise [26]. 

The binary x-constrained flows of the AKNS hierarchy (3.4) are deflned by [15,17,21] 

= -A^i + 9^2, $2,x = + A$2, (3.8a) 
= A*i -r*2, *2,x = -9*1 - A*2, (3.8b) 
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N 



6u j^^Su \bk,+ij '^v<*2,$i>; 



(1) For ko = 0, P = 1, we have 



(3.9) 



By substituting (3.9) into (3.8a) and (3.8b), one gets a x-FDHS [15,17] 
^F^ 9Fi 5Fi 5Fi 

Fi =< A*2, $2 > - < A*i, $1 > + < *2, $1 >< *1, $2 > • 

Under the constraint (3.9) and the FDHS (3.10), the binary t2-constrained flow ob- 
tained from (3.2) with V^'^^ given by (3.5) and its adjoint equation for N distinct real 
number Xj can also be written as a t2-FDHS 

dFo dFo dFo dFo 

g^^, Q^^^ ^l,t2 Q^^. Q^^^ l^-L^; 

F2 =< A^*2, $2 > - < A^*l, $1 > + < *2, $1 >< A^'i, $2 > 
+ < A*2,^l >< *1,*2 > -^(< *2,^2 > - < >) < *2,^1 >< *1,^2 > • 

The Lax representation for the FDHSs (3.10) and (3.11) which can also be deduced 
from the adjoint representation of (3.1) and (3.2) are presented by (2.13) with the entries 
of the Lax matrix M given by [21] 

A{\) = -1 + - E '^^^^Yt''^^' ^ (3-12a) 

2 ^ A An 

BW = tf^- C(A) = f|^. (3.12b) 



A straightforward calculation yields 

^ P 

A\X) + B{X)C{X) ^ P(A) = 1 + Eix^ + JX^tr^' ^^-^^^ 
where Pi, P2JV are independent integrals of motion for the FDHSs 



18 



YUNBO ZENG 



(3.10) and (3.11) 

+ i 52 X ^ \ [(V'lj^lj - '^p2j(f>2j){'^Plk(f>lk - i'2k(f>2k) + ^i'ljhji'2kM, (3.14a) 

^ k^j - 

Pn+j = ^iMij + ^2j(t>2j), 3 = 1, N. (3.14b) 
It is easy to verify that 

N N 

Fi = 5^(A,P, + P^^,.) - ^($]P;•)^ (3.15a) 

N N N N 

^2 = Y.(>^'Pj + '^^jPn+j) - ^^Y.P^) Y.^>^jPj + Pn+,) + lC£Pjf- (3.15b) 

Similarly, it is easy to show that Pi,...,P2n are in involution, (3.10) and (3.11) 
are FDIHSs. The AKNS equation (3.6) is factorized by the a;-FDIHS (3.10) and the 
t2-FDIHS (3.11), namely, if (^'i, ^'s, $i, $2) satisfies the FDIHSs (3.10) and (3.11) si- 
multaneously, then {q, r) given by (3.9) solves the AKNS equation (3.6). In general, the 
factorization of the n-th AKNS equations (3.4) will be presented in the end of section 
3.2. 

(2) For fco = l,/5 = |, (3.8c) yields 

=< *i, $2 >, qx = - < *2, $1 > • (3.16) 
The equations (3.8a), (3.8b) and (3.16) can be written as a a;-FDIHS 

dFi dFi ^ dFi dFi . , ^ 

Fi =< A*2, $2 > - < A*i, $1 > +r < *2, ^1 > < *i, $2 > • 
Under the constraint (3.16) and the FDIHS (3.17), V^^^ becomes 

VrA+ i < ^'i,$2 > - igr J ^ ' 

Then under the constraint (3.16) and the FDIHS (3.17), the binary t2-constrained con- 
sisting of replicas (3.5) and its adjoint system for N distinct real number Aj as well as 
(3.6) can also be written as a t2-FDIHS 

dF2 dF2 dF2 dF2 
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F2 =< A^*2, $2 > - < A^^'i, $1 > +g < A*i, $2 > +r < A*2, *i > 

-\qr{< *2,^2 > - < >) + ^ < *2,^i >< *i,^2 > 

The Lax matrix M for FDIHS (3.17) and (3.19) is given by 

AiX) = -A + i V - (3.20a) 

iJ(A)=,+ if |2f^. C(A) = ,- + if (3.20b) 



A straightforward calculation yields 



J 1 -'A 



A^(A) + 5(A)C(A) ^ P(A) = A^ + Po + Etx^ + (A"^^' ^^"^^^ 

where Po5---5-P2Ar are independent integrals of motion in involution for the FDIHSs 
(3.17) and (3.19) 

Po = ^(< *2, $2 > - < *i, *i >) + qr, 

+ ^ XI — ^[(V'li^lj - V'2i02i)(^lfc01A: - ^2k(t>2k) + 4V'li02i^/'2A:01fc] , 

Pn+3 = ^ ( V'l j '/'I j + V'2j </'2i ) , j = 1 , . . . , AT. 

It is easy to verify that 

TV TV ^ 

Fi = 2j2Pj^ ^2 = 2 J2{X,Pj + P^+,) - -Pi (3.22) 



It is easy to show that Pi, P27V are in involution, (3.17) and (3.18) are FDIHSs and 
commute each other. The AKNS equation (3.6) is factorized by the x- FDIHS (3.17) 
and the t2-FDIHS (3.19), namely, if (*i, *2, g, $1, $2, r) satisfies the FDIHSs (3.17) 
and (3.19) simultaneously, then (g, r) solves the AKNS equation (3.6). 
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3.2 The separation of variables for the AKNS equations . 

(1) For ko = case, we present the Jacobi inversion problem for (3.10) and (3.11) as 
well as for (3.6). With respect to the standard Poisson bracket, ^(A) and B{X) given 
by (3.12) satisfy 

{A{X), Ai^i)} = {B{X), Bi^)} = 0, {A{X), B{y)} = J^^[B{^^) - B{X)]. (3.23) 

In contrast with the B[X) for the constrained KdV flows, the B[X) given by (3.12b) 
has only N — 1 zeros, one has to define the canonical variables Uk,Vk,k = 1, ...,N, in a 
little different way: 

N I , N-1 N 

B{X) ^ = '""""Wxy ^^^^ = n (^-^^)' ^(^) = n(^-^^)' (3-24a) 

j=l ^ ^ ' k=l k=l 

Vk = A{uk), k = l,...,N-l, VN = ^{<^i,^i> - <^2,^2>)- (3.24b) 
The equation (3.24a) yields 

itiv = In < *2, *i > ■ (3.24c) 

Then it is easy to verify that 

{un, B{fj,)} = {vn, A{fj.)} = 0, {vn, un} = 1, (3.25a) 

{uN,A{f,)} = {vM,B{f,)} = B{f,). (3.25b) 

The commutator relations (3.23) and (3.25) guarantee that ui, ...,un, vi, ...,vn satisfy 
the canonical conditions (1.1). Similarly, we define 



VN+j = Pn+j, UN +3 = In— ^, J = 1, A/". (3.26) 

In the same way we can show the following proposition. 

Proposition 5. Assume that Xj,(f)ij,'ilJij G R, i = l,2,j = 1,...,N. Introduce the 
separated variables ui,...,U2n and vi, ...,V2n by (3.24) o-'i^d (3.26). If ui, ...,un-i, are 
single zeros of B{X), then vi,...,V2n and ui,...,U2n are canonically conjugated, i.e., 
they satisfy (1.1). 

It follows from (3.24) that 

Q = e^~, (3.27) 
^2,01, =e-^^, ^ = e«-+^ i = l,...,7V, 

K'{Xj) 1^23 
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or 



^^"^^"""^'^■^ = . = i,...,N. (3.28) 



It is easy to see from (3.13) that 



1 1 ^ 

= -(< $1 > - < *2, $2 >) = X] (3.29) 



Then the separated equations obtained by substituting Uk into (3.13) and using (3.24) 
and the separated equations (3.26) and (3.29) may be integrated to give the generating 
function of the canonical transformation 



N-l . AT AT 



(3.30) 



k=l " i=l 1=1 

The linearizing coordinates are then 

^-1 .u. 



dS 1 1 1 

= 7— = - V / ^=dX- -UN, i = l,...,N, (3.31a) 

aP, 2 (A-A,)/P(A) 2 ^' ^ ^ 



N-l 



'5"+' = 9^ = g / ^137:^7^^''^ + i = 1, .... JV. (3.31b) 

By using (3.15a), the linear flow induced by the FDIHS (3.10) together with the 
equations (3.31) leads to the Jacobi inversion problem 
for the FDIHS (3.10) 

N-l I N 

J2 / 77 ^ , r^ dX - = 7, + (2A, - Y Pk)x, i = I, N, (3.32a) 

^ (A-Aj)V^(^) k=i 

5]/ 7^ -7^^7^7TTd^ + '^N+i = lN+i + '2PN+iX, i = l,...,N. (3.32b) 

t^iJ (A-Ai)V^W 

By using (3.15b), the linear flow induced by the FDIHS (3.11) and the equations 
(3.31) result in the Jacobi inversion problem for the FDIHS (3.11) 



^ rUk 

/ dX — um 

^ I (A-A,)yP(A) 



k=l 
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N N 

= 7i + [2A,2 - Y^^XkPk + XiPk + PN+k) + ^(JZ Pkf]t2, i = 1, N, (3.33a) 

k=l k=l 

N-1 Uk ^ N 

^ 7{ — X \2'r^rT\ ^^ ^ ^^+^ ^ ^^+^ ^ PN+i{4:Xi - Pk)t2, i = i,...,N. 

fe=l (-^ - Xi) VP{'^) k=l 

(3.33b) 

Then, since the AKNS equations (3.6) are factorized by the FDIHSs (3.10) 
and (3.11), combining the equations (3.32) with the equations (3.33) gives rise to the 
Jacobi inversion problem for the AKNS equations (3.6) 



^-1 rUk I 



1 1 '> 



(A-A,)yP(A) 



d\ — itjv 



TV TV N 

= 7i + (2A, - Pk)x + [2A2 - ^^(AfcPfc + XiPk + J + 4 ^fe)']^^, (3.34a) 

fc=l k=l k=l 

^ r>' PN+i ^ 

^ / V\r7mf\^^ + ^^+' " ^^+' ^ 2PTv+ia; + Piv+j (4Ai - ^ Pk)t2, 

(A-Aj) V^(^) k=i 

i = l,...,N. (3.34b) 

If (j)ij,ip2j,q defined by (3.27) and (3.28) can be solved from (3.34) by using the Jacobi 
inversion technique, then (f)2j , ijjij can 

be obtained from the first equation and the last equation in (3.10) by an algebraic 
calculation, respectively. Finally q and r =< ^'i, $2 > 

provides the solution to the AKNS equations (3.6). 

Comparing (2.20) with (3.13), one gets 

JV 

Po = l, Fk = Yli^']-'Pj + {k-l)X';-^P^N+j], k=l,2,..., (3.35) 

where P^j ^ = !> 2, are also integrals of motion for both the FDIHS (3.10) and the tn- 
binary constrained flow. The n-th AKNS equations (3.4) are factorized by the x-FDIHS 
(3.10) and the t^-FDIHS with the Hamiltonian given by 

m=0 li+...+lm+i=n+l 

where li > 1, Im+i ^ 1) Q^m are given by (2.22c). We have the proposition: 
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Proposition 6. The Jacobi inversion problem for the n-th AKNS 
equations (3.4) is 

N-l uu I N 



(A-AO/P(A) 

+2tnY.{-\rar^ J2 \'r'-'Fi^...Fi^^ i = l,...,iV, 



fe=l " V' V " V'V k=l 

n 



2' 

m=0 h + ...+lm+i=n+l 



" 1 ~ ~ 



2' 

m=0 /i+...+Z^+i=n+l 

where h > 1, ...,lm+i > 1? ond Fi^, ■■■Fi^, are given by (3.35). 

(2) For fco = 1 case, with respect to the standard Poisson bracket, A{\) and -B(A) 
given by (3.20) also satisfy the commutator relation (2.31). One define the first N -\- 1 
pair of canonical variables Uk,Vk,k = 1, N + 1, in the following way: 



with 



and 



N N 

R{\) = J] (A - life), K{X) = J] (A - Afe), 
fc=i fe=i 



Vk = 2A{uk), k=l,...,N, (3.37b) 
VN+i = Po = gr - ^(< ^'i, $1 > - < *2, $2 >)■ (3.37c) 
The equation (3.24a) yields 

un+1 — In^- (3.37d) 

Then it is easy to verify that 

{un+i,B{ii)} = {vn+i,A{ii)} = 0, {vn+i,un+i} = 1, 

{uN+i,A{fi)} = 0, {vN+i,B{,i)} = B{,i). (3.38) 

Similarly, we define 

VN+i+j = 2PN+j, j = 1,...,N, (3.39a) 
UM+i+j = ln^, j = l,...,N. (3.39b) 
In the same way we can show the following proposition. 
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Proposition 7. Assume that Xj,4>ij,ipij e R, i = 1,2, j = 1,...,N. Introduce the 
separated variables ui, ...,U2n+i and vi, ...,V2n+i by (3.37) and (3.39). If ui, ...,un, 
are single zeros of B{X), then vi, ...,V2n+i o-nd ui, ...,U2n+i are canonically conjugated, 
i.e., they satisfy (1.1). 

It follows from (3.37) that 

g = e"^+\ (3.40a) 



' 2eUN+i+UN+i+j R{Xj) / 2e^^+i R{Xj) 

^'^^^1 IC{X-) ' ^^^^Y KV^) ' ^ = 1'-'^- 

(3.40b) 

The first separated equations can be found by substituting Uk into (3.21) and using 
(3.37b), the last N + 1 separated equations are given by (3.37c) and (3.39a). They may 
be integrated to give 

S{ui,...,U2N+i) =^{2 \/P{X)dX + 2PN+kUN+i+k) + PqUn+i, (3.41) 
fe=l 

with P{X) given by (3.21). Then the Jacobi inversion problem for the FDIHS (3.17) is 



k=l 

N 



k=l 

N 



/ Pat 
——t^.A + = i=l,...,N. (3.42) 



{x-x,r,/p{xy 

The Jacobi inversion problem for the FDIHS (3.19) is 



N 



k=l 



^^^-^-^^dX^,^ + 2X.t2, ^ = 1,...,N, 



J^ j " _ ^y^^ dA + UN+i+i = iN+i + 2PN+it2, i = l,...,N. (3.43) 



k=i •' 



Finally we have 
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Proposition 8. The Jacobi inversion problem for the AKNS equation (3.6) is 

/life 2 
(A-AOVP(A) ''^ = ^' + ^'^ + ^''-'' ' = ^--'^- 

E / 7^ T^Thttt^^ + ^^+1+^ = ^^+^ + '^PN+it2, i = 1, N. (3.44) 

If (pij, i/j2j, q defined by (3.40) can be solved from (3.44) by using the Jacobi inversion 
technique, then 4>2j,i^ij and r can be obtained from the equations in (3.17) by an 
algebraic calculation, respectively. Finally {q,r) provides the solution to the AKNS 
equations (3.6). 

(3) The above procedure can be applied to all high-order binary constrained flows 
(3.8) and and whole AKNS hierarchy (3.4). 

4. The separation of variables for the Kaup-Nev^ell equations . 
4.1 Binary constrained flows of the Kaup-Newell hierarchy. 

For the Kaup-Newell spectral problem [31] 
take 

"■"^ / „ \2n-2i L \2n-2i-l \ 

K = y<">(„,A)^, V^>i„,X)^j:{j^l^^_,,_, '^^^^^.,,) (4.2) 

i=0 ^ ^ 

where 

ao = 1, a2^~qr, bi ^ -q, Ci = -r, bs = ^{q^r + q^), C3 ^ ^{qr^ - r^), 
and in general a2k+i = &2fc = C2fe = 

fc2k+i\ ^ ^fc2k-i\ a2k = l:d~^{qc2k-i,x + rb2k-i,x), k = l,2,---, (4.3) 

\02fc+l/ \t>2k-lj ^ 

f d — rd~^qd —rd~^rd 
2 V —qd~^qd —d — qd~^rd 
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Then the compatibihty condition of equations (4.1) and (4.2) gives rise to the Kaup- 
Newell hierarchy [31] 

where the Hamiltonian and the Hamiltonian operator J are given by 

7 / ^^ ZT _ 4a2n+2 - rC2n+l - g^>2n+l f C2n+l\ 5-f^2n 



\^2n+l/ 

For 77. = 2 we have 

and the coupled derivative nonlinear Schrodinger (CDNS) equations obtained from the 
equation (4.4) for n = 2 read 

<lt2 = ^(lxx + ^{q'^r)x, rt^ = ~rxx + ^{r'^q)x. (4.6) 

The adjoint Kaup- Newell spectral problem is the equation (2.7) with U given by 
(4.1). We have [26] 

— = =Tr\f "^^"^^ 0iV'2\ dU{u,X) ^ A^i02\ . . 

The binary x-constrained flows of the Kaup-Newell hierarchy (4.4) are defined by 

= -A^^i +5A$2, $2,x = rA$i + A2$2, (4.8a) 
= A^*! - rA*2, *2,a; = -gA*i - A^^a, (4.8b) 



(4.9) 



5m ^ (5m \b2ko+iJ 2V<A*2,$i> 
For fco = 0, we have 

'ci\ ^_(r\ ^l(< A^i,$2> 
,6 J [qj 2V<A^'2,$i> 



By substituting (4.9) into (4.8a) and (4.8b), the first binary x-constrained fiow becomes 
a FDHS 

dFi ^ dFi ^ dFi ^ dFi 

C/Wi OW2 0^1 0^2 
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with the Hamiltonian 

Fi =< A^*2,^2 > - < A^*i,$i > -7: < A*2,^i >< A*i,$2 > • 

Under the constraint (4.9) and the FDHS (4.10), the binary t2-constrained flow ob- 
tained from (4.5) and its adjoint equation for N distinct reral numbers Xj can also be 
written as a FDHS 

dF2 ^ dF2 ^ dF2 ^ dF2 , , 

with the Hamiltonian 

F2 = - < A'^*2, $2 > + < A"^*!, $1 > +^ < A*2, $1 >< A^*i, $2 > 

+^ < A^*2, $1 >< A*i, $2 > < ^*2, 3>i >^< A*i, $2 >^ 

+ J(< A^*2, $2 > - < A^^i, $1 >) < A*2, $1 >< A^'i, $2 > • 
8 

The Lax representation for the FDHSs (4.10) and (4.11) are presented by (2.13) with 
the entries of the Lax matrix M given by 

B(X) = lxt^0^' C(A) = iAf^|^. (4.12b) 

j=l 3 j=l j 

A straightforward calculation yields 

^ p ;^4p2 
^^(A) + 5(A)C(A) ^ P(A) = 1 + E[a^ + (A^^^^' ^^"^^^ 

where Pj, j = 1, ...,2N, are 2A/" independent integrals of motion for the FDHSs (4.10) 
and (4.11) 

11 1 

Pj = -2^|(V'2i02i - V'lj'/'lj) + g < A*2, *1 > ^ji^ijhj + g < A*i, *2 > >^ji^2j(t>lj 
1 . 1 

+ 8 T2^rT2 [''^?'^fe(V'lj'?!'lj-V'2j02jO(V'lfe01fe-V'2fc02fc 



28 



YUNBO ZENG 



j = l,...,N (4.14a) 

Pn+j = \ {ipij (1)1 j + (t>2j), j = l,...,N. (4. 14b) 
It is easy to varify that 

N N N 

Fi = -2 5^ P„ F2 = 2 5^(Af + A|P^^,.) - -(5^ P,f, (4.15a) 



2 



N 

< *2,^2 > + < >= 4^Piv+j. (4.15b) 



By inserting A = 0, (4.13) leads to 



l + \{< *2,$2 > - < >) = /P(0) 



H-5][-P,A-2 + p2^^.]- (4.16) 



With respect to the standard Poisson bracket it is found that 

{^(A),A(/.)} = {S(A),i?(/x)}, (4.17a) 

{A{\),B{^)} = 2(A2-^2) [/^-gM - ^^(^)]- (4-17b) 

Then {^^(A) + 5(A)C(A), ^2(//) + P(/x)C(/x)} = imphes that Pj,j = 1, ...,2iV, are 
in involution. The CDNS equations (4.6) are factorized by the x-FDIHS (4.10) and 
the t2-FDIHS (4.11), namely, if (^'i, ^fa, $i, $2) satisfies the FDIHSs (4.10) and (4.11) 
simultaneously, then (g, r) given by (4.9) solves the CDNS equations (4.6). The factor- 
ization of the n-th Kaup-Newell ewuations (4.4) will be presented in the end of section 
4.2. 

4.2 The separation of variables for the Kaup-Newell equations. 

Since the commutator relations (4.17) arc quite different from (2.31) and (3.23), we 
have to modify a little bit of the method presented in sections 2 and 3. Let us denote 
A = A^, Aj = A|. The entries of the Lax matrix M given by (4.12) can be rewritten as 

A(A) = l + ^(< *2,$2 > - < >) + ^AA(A), S(A) = ^VaP(A), (4.18a) 
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where 



N 



N 



Xjll^2j(t)lj 



It is easy to see that 



{A{X),A{Ji)} = {B{X),B{Jl)} = 0, 
{ACx),B{fl)} 



^ ^B{fl)-B{X)]. 



X — jji 



(4.18b) 

(4.19a) 
(4.19b) 



It follows from (4.16) and (4.18a) that 



A{X) 



N 



(4.19c) 



The commutator relations (4.19) and the generating function of integrals of motion 
(4.13) enable us to introduce ui, ...,UAr in the following way: 



TV 



B{X) = E 



K{xy 



(4.20a) 



with 



N-l 



N 



RiX)= l[{X-Uk), KiX) = l[{X-Xk), 
fe=i fe=i 

and vi, vn by ^(A): 

Vk = A{uk),k= 1, ...,N - 1, 
Wat = - < *2,$2 > ■ 



The eq. (4.20a) yields 
Similarly we define 



itAT = In < A*2,^'l > ■ 



VN+j = 2PN+j, j = 1, N, 

UN+j=i.n—, j = l,...,N. 



(4.20b) 
(4.20c) 

(4.20d) 

(4.21a) 
(4.21b) 



Then we have 
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Proposition 9. Assume that Xj,4>ij,ipij e R, i = 1,2, j = 1,...,N. Introduce the 
separated variables ui,...,U2n and vi, ...,V2n by (4-20) and (4.-21). If ui, ...,un-i, are 
single zeros of B{\), then v\,...,V2n and ui,...,U2n are canonically conjugated, i.e., 
they satisfy (1.1). 

It follows from (4.9), (4.20a), (4.20d) and (4.21b) that 




(4.22a) 



By substituting Uk into (4.13) and using (4.16) and (4.19c), one gets the first N — 1 
separated equations 



Vk = A{uk) = —[Jp{uk)-^/P{0)], k = 1,...,7V- 1, (4.23a) 
where -P(O) are given by (4.16) and 



P(A) = 1 + V[.^+ J ]■ 
A-Af (A-Af)2^ 

It follows from (4.15b), (4.16) and (4.20c) that 

N 

VN ^2-2^/Pi^-2J2PN+^■ (4.23b) 

i=l 

The separated equations (4.23) and (4.21a) may be integrated to give the generating 
function of the canonical transformation 

Siuu...,U2N)^ J2U ^\/Pi^)dX-2y/P(^ln\uk\] 
k=i A 

N N 

+ (2 - 2/P(0) - 2 ^ PN+i)uN + 2J2 PN+iUN+i. (4.24) 

1=1 i=l 

The Jacobi inversion problem for the FDIHS (4.10) is 

/-"fc 1 ~ 1 1 
[ / dX H /nlttfcl] H um = 7i — 2x, 
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fe=i A(A- Af)2-^P(A) V^(0) 



PN+i 



-{^== + l)«iv + «iv+i = iN+i, i = 1, iV. (4.25) 
The Jacobi inversion problem for the FDIHS (4.11) is 



7V-1 



/-Mfc I ~ 1 

Et/ 7^dX+ ln\uk\] + 



Mat 



AT 



7i + (2A2-^Pfe)t2, 



fc=i 

fc=i A(A- Af)2^P(A) V^(0) 
-( + l)^jv + UN+i = 7iv+» + 2XtPN+it2, i = l,...,N. (4.26) 

Finally, since the CDNS equations (4.6) are factorized by the FDIHS (4.10) and 
(4.11), combining the equation (4.25) with the equation (4.26) gives rise to the Jacobi 
inversion problem for the CDNS equations (4.6) 



y^[/ — ; dXA 1 Inlukl] H 1 



TV 



7, - 2a; + (2A2 - J] P,)t2, i = 1, iV, (4.27a) 



fe=i 

-^-1 cUk \4p p 



fc=W A(A-Af)V^'(A) V^^°) 



-{^== + 1)UN + WAT+i = 7iv+i + 2A^Piv+it2, i = 1, AT. (4.27b) 

If (t)ij,ip2j,q defined by (4.22) can be solved from (4.36) 
by using the Jacobi inversion technique, then (f)2j , V'l j can 

be obtained from the first equation and the last equation in (4.10), respectively. 
Finally q and r = — < A^'i, $2 > provides the solution to the CDNS equations (4.6). 
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In general, the above procedure can be applied to the whole Kaup-Newell hierarchy 
(4.4). Set 

oo 

A2(A) + 5(A)C(A) = ^FfeA-2^ (4.28a) 

k=0 

where F^, A; = 1, 2, are also integrals of motion for both the x-FDHSs (4.10) and the 
tn-binary constrained flows (2.16). Comparing (4.28a) with (4.13), one gets 

N 

Fo = l, Fk^J2^Xf-^Pj + {k-l)XfP^+j], A; = 1,2,.... (4.28b) 

By employing the method in [28,29], the in-FDIHS obtained from the t^-constrained 
flow is of the form 

dFn ^ dFn , dFn , dFn . , , 
with the Hamiltonian 

n— 1 ^ 

^n = 2j:(-^)-^ Y: F,,...Fi^^„ (4.29b) 

m=0 h+...+lm+i=n 

where /i > l,...,lrn+i ^ 1; and am are given by (2.22). Since the n-th Kaup-Newell 
equations (4.4) is factorized by the x-FDIHS (4.10) and the t„-FDIHS (4.29). We have 
the following proposition. 

Proposition 10. The Jacobi inversion problem for the n-th Kaup-Newell equations 
(4-4) ^5 given by 



[ / dX H /nlttfel] H um = 7i — 2x 



+2*»E(-^)"'""' E Af"+'-'F,,...F,„, i=l,...,N, (4.30a) 

m=0 li + ...+lm+i=n 

El XfPN+i ~ PN+i , I II / Pn+i , , 
[/ — 7^^^^ T^T^^^I^felJ - ( ^—- + l)UN + UN+i=lN+i 
fe=W A(A-Af)yp(A) vno) vno) 



n-1 



+2tnJ2{~ram J2 {lm+i-l)Xf'^^'PN+iFi,...Fi^, z = 1, TV, (4.30b) 



2 

m=0 li + ...+lm+i=n 
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where li > l,...Jrn+i > l, Fi^, ...Fi^, are given by (4-28b). 

For example, the third equations in the Kaup-NeweU hierarchy with n = 3 are of the 
form 

13 13 

Qts = -^Qxxx - g(9^''^ + ^QfQx)x, ^3 = -^r^a^x - gi'^^Q^ - '^(lf"rx)x- (4.31) 

The Kaup-NeweU equations (4.31) can be factorized by the x-FDIHS (4.10) and t^- 
FDIHS with the Hamiltonian F3 defined by 

TV N N ^ N 

F, = J2{2X'^Pj + 4A5P2 ^,.) - E(A|P, + A,^P2 ^,.)] E + 4 (E P^)'- (4-32) 
The Jacobi inversion problem for the equations (4.31) is given by 

/-Mfe 1—1 1 

iV AT 

= 7i - 2a; + [2Xt Y.^X]Pj + A^P,. + A^^P^+^O + -(^P,f ]t3, i = 1, 

k=iJ A(A-A2)yp(A) V^(0) V^(0) 

+MAr+i = 7Ar+i + [4AfPiv+i - A^Pjv+j- E ^:?]^3, i = 1, AT. 

i=i 

In general, the method can be applied to all high-order binary constrained flows (4.8) 
and whole KN hierarchy (4.4) in the exactly same way. 

4. Concluding remarks. 

For high-order binary constrained flows, the method in [1-6] allows us to directly 
introduce N + kQ pairs of canonical separated variables and N + ko separated equations 
via the Lax matrices and generating function of integrals of motion. In this paper 
we propose a new method for determining additional pairs of canonical separated 
variables and N additional separated equations for high-order binary constrained flows 
by directly using AT" additional integrals of motion. This method is completely different 
from that proposed in [23,24] and can be applied to all high-order binary constrained 
flows and other soliton hierarchies admitting 2x2 Lax pairs. 
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